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ABSTRACT 


The basic equations of quantum scattering are translated into the 
Wigner representation. This puts quantum mechanics in the form of a stoch- 
astic process in phase space. Instead of complex valued wave functions and 
transition matrices, one now works with real valued probability distributions 
and source functions — objects more responsive to physical Intuition. Aside 
from v^ritlng out certain necessary basic expressions, the main purpose of 
this paper is to develop and stress the interpretive picture associated with 
this representation and to derive results used in applications published 
elsewhere. The quasi-class ical guise assumed by the formalism lends itself 
particularly to approximations of complex multi-particle scattering problems. 
We hope to be laying the foundation for a systematic application of statisti- 
cal approximations to such problems. Tlac form of the integral equation for 
scattering as well as its multiple scattering expansion in this representa- 
tion are derived. Since this formalism remains unchanged upon talcing the 
classical limit, these results also constitute a general treatment of classi- 
cal multi-particle colllsiun theory. Quantum corrections to classical pro- 
ppgators are briefly discussed. The basic approximation used in the Monte- 
Carlo method is derived in a fashion which allows for future refinement and 
which includes bound state production. The close connection which must exist 
between inclusive production of a bound state and of its constituents is 
brought out in an especially graphic way by this formalism. In particular 
one can see how comparisons between such cross-sections yield direct physical 


Insight into relevant production tncchnnistns. Finally, as a simple illustra- 
tion of some of the forroallf-a we treat scattering hy a bound two body system. 
Simple expressions far single and double scattering contributions to total 
and differential cross-sections as well as for all uecc^ssary shadow correc- 
tions thereto, are obtained. These are compared to previous results of 


Glauber and Goldberger 


I, INTRODUCTtON 


This paper developes tlu-j uluments of whut may be called the Wigner 
representation of quantum scattering. It is well kno^m that quantum mech- 
anics can bo formulated entirely in terms of density opezutors and, linear 
operators (often called superoperators in this context) which act upon the 

n 2 9 1 

densities. J The Wigner representation of densities puts them in the 

form of real functions of the coordinates and momenta of the system's part- 
icles. In other words, density operator matrix elements are functions 
defined in the classical phase space of the system. Schroc* linger 's equa- 
tion, correspondingly transformed, is a linear equation for the density's 
time dependence. It loolcs like the equation governing a Markoffian 
Stochastic process in phase space. The integral form of Schroedinger 's 
equation, the Lippmann-Schwinger equation in the usual Hilbert space repre- 
sentation of the theory is transformed into an integral form of the equation 
of continuity. The equation of continuity relates densities, currents and 
sources. Thus, upon translation into the Wigner representation, the funda- 
mental equations of quantum scattering appear in terms of these physically 
graphic and transparent objects. 

The main purpose of this paper is to develop and stress this pic- 
ture and, in particular, to indicate hov; its graphic, intuitively appealing 
nature may be exp.latj.ted Ir* the formulation of approximations to complex 
multi-particle scattetlag problems. 



The I'oCTialiBm Mn.oini>af nen > -jUi las., * ■:’ a' ‘i quantun Hcattorlnp. 

In fact the Wigner veprcbentatlon iu c-.ij)C' f.illv conr.ti*', '.ipd to pc. over 'mooth" 
ly to the classical limit. In taking thi;; It nit none >l the cqoation.% change 
tUelr form nor do any of the symbols ‘n their change in iuterpv'tta- 

tion or role. Schroedinger 's equation goes o’'cr into the Llonvllle equation 
and its integral forri niinilar fy cn:ei‘f;ca 'is the >■. r'.sslecil let'/prol cquat en t'f 

n,) I 

motion obtained here in its most general foiiuy* 

Quantum mechanics assumes a classier.’. ',.’i ■ rn thi' fo^'m.illsmj in 
this lies the source of its intuitive ^ .vc'-r. MevevtiH'less . since we have 
here merely another representatiim of oidi-v'iv quantum i. et J- .ni»! , nil of the 
latter's wave and interference properties must bs >it<lden vithin. Wmt has 
happened is that this translation of quantcjm mrcuaiiict} iiecer;.b‘ari]y gives 
birth to a non-classical stochastic process. Densititur are not positive 
definite and in addition exhibit long range ''scillatory bchavlor--iust that 
needed to produce Interference!. Ihus although the total structure of the 
formalism developed here and the various roles played V'y its dcT’ents (densi- 
ties, sources, etc.) are i&nnorphlo to cl.asslccl theory, tlv particular 
functions needed to represent these eli'nv.jncs will differ conJidcjrably in some 
respects from those of classical stochastic theory. 

These qualitative remarks point to the ciroumstavicps in which it 
may or may not be useful to en.ploy this representation. The property being 
observed should not depend, critically on Hich or *er Interfereucu effects or^ 
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concomitantly on the existence o£ certain long range order or correlations. 

There also should be a significant advantage in being able to visualize the 
process quasi-class icallv. By this we monn tliat during the time in which 
the reaction takes place, at least some of Its main participants can be use- 
fully pictured as having simultaneously defined positions and momenta while 
moving along trajectories in classical phase space. 

Generally these properties occur in conjunction and typify systems 
which have been difficult to treat by other methods. Thus one may hope tlj;i-c 
the approximations engendered by this representation will be complimentary to 
those previously established. 

.•4‘ 

An example satisfying these conditions occurs in medium and high 
energy collisions between various projectiles and nuclei. Elastic scattering, 
most of it diffractive, is treatable by the Glauber approximation, the multiple 
scattering expansion, the optical model or, various combinations of these. This 
is a highly coherent process depending critically on high order interference 
effects. Inelastic scattering, in which the target is left in one or possibly 
a few well defined final states, may also be treated by such methods. This 
leaves about half of the possible reactions unaccounted for — the non-elastic 
collisions — in which many nuclear species may be produced in assorted multi- 
plicities and momenta. Since these final states are so complex, one generally 
observes averages over them such as inclusive cross-sections and multiplici- 
ties. This washes out most high order interference effects. Since this for- 
malism is v/ritten directly in terms of the density matrix, such averaging 
may be automatically done with utmost case and elegance. In contrast, when 
calculating in the usual Hilbert space with pure final states, one must en- 
umerate them, square and approximately sum over unobserved final states (e.g. 


using closure). Such procedures are finessed by this fonnnllera. Finally, 
perhaps the most important attribute this formalism brings to such e, problem 
is its quasi-classical guise. These collisions involve mruiy particles, high 
orders of multiple scattering, possibly collective motion, and large oergy 
and momentum transfer. It is vor; dlffictilt to ftoe sv.ch phenomena ill rarms 
of waves in multidimenaicnal position or ntum space, ilu' minds cy i;i 
irresistably drawn to hydrodynaiisic, thermodynamic or transput i. ; iicoi'y type- 
pictures in which joint average distribution functions in puiutlon and mom- 
entum play a central role. In contrast to previous ad hoc application of 
such concepts to such problems of scatt-^^ring, they occur here .as natural and 
systematic approximation procedures of ordinary quantum mechanics. 

Another field which may be mentioned is chemical reactions. Hare 
again complexity Is most often the rule, statistical averaging the natural 
ally, the Important steric properties of the compounds Is easily represent- 
able and, an added positive factor, short wavelengths often mahe the problem 
semi-classical as well as quasi-classlcal. The Wigner representation is 
especially suited to the semi-classical limit. 

fill 

There are other potential .areas of application*' ‘ ^ whicli will not 
be discussed since those already mentioned should serve to illustrate practi- 
cal reasons and requisite criteria for use of the reprwtnieticn. However, 

it should be finally noted that the Wigner representation has bad a long 

fl2l 

history of use in, especially, transport theory. l*Jlift is being developed 

here is a version suited to collision phenomena cn complex yet microscopic 
sys terns . 

A synopsis of the remainder of this article iollmiJ.a. 
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Section II: Elementary definitions, exemplos and theorems connected 

to the Wigner representation are given. A Bra-Ket notation is introduced to 
designate elements of the vector space of density operators. The mapping of 
operators from their usual representation into the nev one is given* 
Schroedinger ' a equation is then transformed into the new representation. 

Section III: The serai-classical limit is very "briefly discussed. 

The propagator for finite time translations on densities (not pure states) is 
obtained as a perturbation series in terms of its classical, not free, value. 
Thus, just as one thinks of waves propagating freely between successive scat- 
terings by a potential in the usual repr-. ;Tentation of perturbation theory, 
here a particle moves on its classical trajectory between successive quantum 
jumps. These ideas lead to practical formulae for computing quantum correc- 
tions to semi-classical processes but this branch of the subject will not be 
developed in this paper. 

Section IV; The most basic concepts used in scattering are intro- 
duced in the context of potential scattering. The time independent integral 
equation for scattering is derived. The effect of the potential on the in- 
coming si.i?.tionary flow of particles from the accelerator is expressed by a 
source distribution function in phase space. This gives the net rate of pro- 
duction of particles being produced by the potential with a certain position 
and momentum according to the heuristic, quasi-class ical interpretation of 
the symbols in this formalism. Source functions are in many ways analogous 
to transition matrices but here occupy an even more central role in the theory. 
The relation between the solution to the integral equation and the obser/ed 
scattering cross-section is established by noting that the spatial integral 



c 


of the local production rate of particloo of a certain jnojnontum io the total 
production rate which in turn io the incident flux timeo the differential 
crooo-section. 

Section V: Tte reoulto of the previouo flection are extended to the 

multiparticle problem where, it io expected, the moot effective uoo of the 
formalism io to be found. Formulae continue to be analogouo to ctandard re- 
QUlea of scattering theory in multiparticle Hilbert space. 

Section VI; The multiple scattering expansion is derived. The re- 
sult is puite fllmilor to that of VJatson. Here the role of the fully off 
shell transition matrix is token by a 'jump operator'. 

Section VII ; The general properties of the jump operator for two 
particles interacting via a phenomenological potential are fully as mysteri- 
ous as those of its counterpart, the transition matrix. Much energy has 
been expended but little physical intuition can illuminate the off-shell prop~ 
erties of the latter complex function. Two limiting forma of the jump oper- 
ator are examined, the classical and the dilute. The relation between a 
particle moving on a classical trajectory and the concept of sources previously 
developed is explained. In the dilute limit, in which all particles are very 
far from each other compared to all other length parameters of the problem, we 
obtain the basis of the Monte Carlo method. However, having obtained this 
not as an ad hoc proced\ire but as a well defined quantum mechanical approxi- 
mation, one sees immediately how its range of application may be properly ex- 

f5l 

tended to include bound final states and in addition, how systematic 
improvements in the approximation may be made for less dilute systems. 

Section VIII: The multiple scattering formalism is developed 

further and then applied to find expressions for various simple but important 
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production croso-sectiono. In addition, the close relation tetwoen inclusive 
production of a bound state and of its unbound constituents is discussed. It 
is pointed out how comparison between these mcaoxweraents can provide informa- 
tion on production mechanisms. 

Section IX: As a final illustration of the formalism we look in 

some detail at the simplest problem to which it might usefully be applied; 
scattering of an elementary projectile by a two body bound state in the di- 
lute limit. The lowest order term in the multiple scattering expansion 

gives immediately both the differential cross-section formula first derived 

F61 

by Goldberger for quasi-free scattering by a bound particle^ ^ as well as a 
corresponding expression for the total cross-section. The second order term 
gives an obvious shadow correction to Goldberger *s formula (the effects of 
which may be seen quite dramatically in reducing the backward peak in 
proton-deuteron elastic scattering at intermediate energies), a shadow 

r?! 

correction to the total cross-section which is compared to that of Glauber'' 

and, the expected double scattering contribution to the differential cross- 

section. All these results are straight-forward to obtain and transparent 

in physical meaning. They are written in terms of observed cross-sections. 

They are extendable in a straight forward fashion to more complex systems 

f81 

and also to non-dilute systems. ■ ^ 


II, ELEMENTS OE THE WIGNEK lUiVlUJSF.NTATION 


We begin with o jrc?iumG ot previous reau3» 6 concerning the Wlgner 
repreaen Cation. Details oC certain atandarcl dcrivationn which mav be found 
in the earlier literature* .,'co omitted. Ini tl; V', di -.^usoion will he confined 
to the case of one apinless particle. Tlic extension to more than one spinleso 
particle is straight-fomord and v»ill be used in later sections. Throughout 
the paper we ignore spin and Indtstinguish.ability of particles. 

Let C5^be any opetator o.i single particle Hilbert space. Its 
VJigner representative is defined to bo^^^ 


• <21 - » J. I ‘ > 

I d ^ e ' ^ ^ 


CO 


Wc will givG Qome glmplo examplen. Let 2^ end denote the ordinary position 
and momentum operators and let ) denote some functions of these opor- 

ators ; then 


/ = (9'C30 


and 

( (9'CP))^ (£, = <y-c 


C^> 
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Wlien depends on pvoducts of conjugoto oporacors* then their ordering will 
hove duo c£l«ct» Another instructive example is the density operator of a 
normalized Gaussian wave packet state in Hilbert space* The Wigner repre- 
sentative of this operator is 

IT a, L(. 2 

I ’ 

- 1 a - F") ^ 

where and are averages and /S 
For any operator 

(ITT)-’ i Z. , -ff. ) , (4) 

|<9| = j dx <Z \ (»l 

f dx &rr)-^ (XjCx,^K 

Xn particular, if iP- is a density operator jCf , then the normalization con- 
dition is 

^ I - 1 ' 

It can be seen that appears to be a sort of Joint probability 

distribution function in phase space, the volume element of which is given 
by 

c/ ^ tiTi) - d x 

(since "ti “ i ) . We find it convenient to denote a point in phase space by 
a single 33onbol 

if = ! J 



iO 


and correspondingly q volume cletncut no 

dcf « 6 1C ? 

It Is straight fon^ard to slicA/ In gencrol tl^at for any two operotora cP^ 

and 

{(9-C9^'\ = f ‘‘H’ “f > 

In particular, the expectation value of any in the .itatu dciicribed by 
is given by 

'riiia is consistent with a probability t'ciisity interpretation of pij . Ivliat 
is inconsistent v;itli this Interpretation hov/evor, at least in Hio 

ordinary sense, is the faet that altiinugh is necessarily real, it la not 
necessarily positive everywhere. Thus 

<< M- ) 

carjK.t 'cfer to the ptoboblllLy of a roali/a’vle mearujrement (measuring whether 
the particle is In tor arbitrary regions, pliaso space. This 

is consistent with the 'luaiUntm ricdianical fact that to every rfiglon of phaec 
space there does not neci.oaarily correspond a physicolly rcall^sable mensure- 
ment. In particular, A mi.dit violate the uncertainty principle limit. 

As wc shall sec, this foc.i iloes in t appear to affect the. heuristic value in 
thinking of as a joant probability distribution function in the quantum 
8cattei.li',g formalism to be developed. 

liquation 1 maps Illlhcrt space operators onto pbasc space functln'as. 

Tlie inverse mapping may be performed with the help of tlie operator 

R tx , - 5 d (?. e ’ \ ^ ti 

which has Wlgner rcprei;entativo 

(-i , L ;£') = (iti)’ S(2-z'> SiE-E'^ 

(lO) 


( 9 ) 
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W) )tj <" ) 

Thus maps Into the Dirac delta density in phase space. Although it is 
not a true quantum mechanical density operator, since it does not describe a 
physically realizable state of the system, it will be a convenient heuristic 
device to speak of it as a density. 



We now note some useful properties of 


J = 1 , 

) 

\ cifM =■ 1 , 

Cl3 ) 

= 1 (9- R tif) 1 ^ 

) 

1 fU*f) 1 = SLUf, Cf') . 

0 5- ) 


The above equations imply 

c9- = ^ \ <9- I 


U*- ) 


In fact it is easy **,o see that this integral has the same Wlgner representa- 
tive as C?: It therefore must equal C^as long as distinct Hilbert space 
operators map into necessarily distinct phase space functions. Tlie latter 
statement is true because the Wigner representation is accomplished via a 
Fourier transformation on the matrix elements of (9r 

Another bit of notation is now introduced which serves, ai ng other 
things, to enhance the analogy between the usual quantum scattering formal- 
ism and the representation of it to be developed here, is considered 

[21 

as an element of a linear vector space and written in bra-ket notation as 


\( 9 > . 


22 


) (5- > -a ( 9 ^ Cl7) 



•55 1 1 

= J C9, i*f) (Sh, (cp) 


We also write 

1 R(.q>)> = 

1 (f> , 

(n ) 

from which one gets 

<t(>l ’ 


C7.r> ) 

< <f 1 s>'> •= 

"5 C tp^ S’ ' ) 

CzO 

j d S’ 1 M’XS’i = 1 



Quantum operators appear in this formalism in two distinct ways; 
as kets, of the type just introduced, or as operators on kets. Generally, 
density operators are mapped into kets. Now let (5^ and ^ be any Hilbert 
space operators (yO generally being a density operator) , we write 

I yo "> ® 1 > ! ' 

\ p'y ^ \ /o <9-> ■ fz3 ) 

f 31 

Previous authors have shown, *• •' in effect, that the matrix elements of these 
operators are given by 

< If l I If ' > = ) ( e-MJ ^ < ‘f 1 M>’> 

< cp I tp'> - c9ij (S’) (e^-fn 1 <‘f • S’' 

^ ^ (34) 
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where A Is the Poisson Bracket operator , 

f Ag ^ a)- ■) 

Alternatively we note 

= I P.(v) (^fUun \ . (z^) 

r 

/, 

f 

Similarly it is easy to see that ■ 

<«-pl (5»t I cf '> = \ Kc<v*) <3- I (^3-j ) 

Thus they are related by transposition* 

< <f \ \ Lf'y ^ < cf U (^i. I 

j 

= <r I <9 -r t tf U S’) ; 

We are now prepared to consider the form which dynamics Cakes in 
the Wigncr representation, Schroedingcr *s equation is 

- ^ H') 

as written in the usual Hilbert space# In terms of Wigner representatives 
this becomes 

^/hiz \ py = '■i i H > I 

|| 

^ C 2-0 ) i 

I ; 

It is easy to see from the preceeding that the hermeticity of W implies | 

(.30 ll 


Li 


from which are abtainecl 


D = . C.S 2. 

In addition we note 

D 


The antisymmetry of the time evoli .ion operator is sufficient to prove con- 
servation of probability and time reversal invariance. 

By choosing to think of <C. as a nrobabllity, this repre- 
sentation of Schroedlnger *s equation appears to describe the time evolution 
of a 'stochastic process', . Schroedinger's equation becomes its 

differential equation. The procedure leading from this equation to the mas- 
ter equation and thence to regions of statistical physics has been under 
investigation for a long time. This paper is not an attempt to contribute 
to such lines of research. Rather, we wish to show that the analogy to 
statistical physics implied by the probability density Interpretation of 
^Cj’l^^raay be carried through to all the basic equations of quantum scatter- 
ing theory. In fact quianUum scattering theory falls most naturally into the 
mold of a stochastic process ‘kiscribablu iu classical phase space. 
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III, QUANTUl-t CORRECTIONS TO CLASSICAL DYNAMICS 
The dynamical equation can be re-expressed as 

where Hi<f) Is just the classical Hamiltonian. The first term In the expan- 
sion of the sine gives the classical liouvllle equation 

V?>i: = (3 A") 

Let US call operator which generates classical time evolution; 

It Is clear that the solution to the classical problem is given formally by 

^ ( 37 ) 

where Is the function which gives the coordinates of the phase point of 

the classical trajectory at time t such that Cfa - <fJ 

Let the operator which corrects the classical propagator, 

e - Mtr) 

It satisfies 

%^Ate)= e. C D-x>c') e " 

= A Ct» -■!)<.■> e“ C3-?) 

Its perturbative solution with initial condition t is 

A lir) - I + dt=‘ <2^‘*'CD-Pci e"'’"* 

C 4 o) 


+ -- - 
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The first quantum correction to the matrix element of the time evolution 
operator is thus given by 

= <•<=' < I 0-Dc» V 

using equation (37) . 

The interpretation oE this equation is that a particle may travel 

I 

from to in a time til via an infinite variety of paths in phase space* 
Each path is itself made up of some nti..iber of classical path segments, the 
first beginning at and the last ending at . The particle moves class- 
ically along the first segment, performs a quantum Jump through phase space 
to the second when it reaches the end of the first and, proceeds in like 
manner until the end. Each quantum jump, say that going from point to 
, occurs with a probability D- ‘Oc i . This probability 

may be negative. In fact, due to the antisymmetry of both D and Do , the 
time-reversed jump always occurs with a sign opposite that of the time- 
direct Jump, 

It is interesting at this point to exhibit in more detail the 
jump probability operator's matrix elements for potential scattering. If 

^ p*/ 2 iv^ + ^ ^ ^ contri- 

bute to D- Pc . We have 

^ hy Jy' ^ - f'jt' 

' I I ^ ^ if > < ^ i V » Jf - O'l I 

^ 5 CTc -X') \/(x~i y ) 
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so that 

D-Pc Ilf'?' 

= SC2-Z0[P^ 

' 1 ( J> - VtZ-i ?>) %H. V/tx) ^/i>^^lf-«'>)\ 

ei3) I; 

This shov;s that a particle may experience a finite jump In momentum but not I 

in position, which Is what is expected for a particle subject to a random 
force. One can show that the average of this random force Is just the 
classical force. 



IV. CURRENTS, SOURCES AND CROSS-SECTION.S IN POTENTIAL SCATTERING 

In this section wc wish to develop the essential Ideas of this 
paper in the simplified context of potential scattering. These results will 
subsequently be extended to the general multiparticle, multi chnmiel case. 

Tlie plan here is first to derive the integral equation for scattering — the 
Lippmann~Schwlnger equation in the Wigner representation. Next to relate 
the solution of this equation to the cross-section. Finally to interpret 
the symbols appearing in these results in a way which brings out clearly the 
analogy to statistical physics. 

Let 

be the time evolution operator in the absence of interaction. Explicitly, 
using the expansion in terms of Poisson Bracket operators, it is 

< tf n>o * - i!! • H 5' 3 

where 

Let represent any density of freely moving particles, 

V7>t: 

and consider a density satisfying 

I ^ \^Lt^y 

C^7) 
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then clearly 


and 


JLUrx. 

~ <o 


( I <.^r) *> - \ p ^ ^ « O . 


m) 

w) 


TIjub is that solutior to the potential scattering problem with 

boundary condition that it approach a specific free density as 

««• X 30 • 

We next specialize to the case of stationary flow. If ^ It) is 
chosen to be time independent then t) must also be stationary and 

the integral equation becomes 

I ^ l/^> + [ i } (D-D.*) 

Tlie upper limit on the integral, being arbitary, is chosen for convenience 


at f:- D . Now since 
^ X>«^r 


1 Z, ^ \ Tt t, 42 > 


then 


-T>ot 


(St) 

(5£) 


<: -i. \ e.' = < X *r *=, iP I 

Therefore the integral equation can be rewritten as 

< ’ C2£, i/’> 

+-f d* < ij +jrt, I .C-5'i) 

-CM ' 

As wc shall see shortly, this equation has a very simple interpretation. We 
note in passing that 




t e 


- "D o t 


— dO 




- CT>o - 7 - ^ t 




= - ( p„ 


(s4) 


(£■^1 


In the limit , TIujb wg enn olso write 

Xp'^'y = \f>y - ct>o-7) ( D-T).,") 

To obtain on expression for the cross-section \m first consider the 
standard incoming beam of one particle per unit voJume. Thus v;e set 

\ f>'> 

llie state represented by this density provides an entering flux - . One 

denotes by / the corresponding solution to the stationary scattering 

problem. Tlien 

Tgiz, e‘*>> 

is the current density of particles of momentum ^ sot up by the interaction 


^ ) 


) 


and 



where /S is a surface totally enclosing the region of interactionj la the 
total outgoing flux of particles within thu interval about ^ . Let 

ert^p, be the probability that an Incident particle of momentum 

^ ^ be scattered by the potential into the corresponding Interval. Thus 

CT( -p, (.tn-)-’ (ss) 

Applying Gauss’s theorem we get 

CT(^, - {A-vy^S S.r« 2, or?) 

where the integral extends over all space and the source function is 

defined by 


C Tf , ^ ±G (.Tf 

= “ <H) ^ i i>o i 

^ I D-Do \ (to) 
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The equation of continuity wlion applied to a stationary i -tributlon 
shows that the divergence of the current> the source function as dctiu<’.d here, 
does indeed give the local rate of production of particles within the element 
«• cf , Our equation for the cross-section therefore 

merely states that the total rate f proauction of particles in interval 
is the integral over all space of the local production rate. It is easy to 
gee that the usual differential cross-section is given by 

<r(|L, ^e) = I (6 1 ) 

The equation for the sour.,e function in terms of the density func- 
tion can be rewritten as 

Sg. c ^ <<(>\ P-'Do \ ' I 0 ^^ ) 

Now is the total probxibility per unit time that a particle 

at will jump to , while <T I is the contribution to that 

probability due to simplr free particle streaming. Thus by the above equation, 
since if * ^if'l .*s the number of particles in steady state in \ 

S^tf) is the rate at which particles are jumping into due to the 

potential. Tills is equal and opposite to the rate at which particles are 
jumping into tfef due to streaming. The net rate at wfliich particles are jump- 
ing Into Jcf is zero since we are in steady state i.e. 

i d cf' <q>! r> f cp' yCif'/ ^ a, 

The Wigner representation of the Llppmann-Schwrlnger equation for 
the time independent wave function of the scattering problem can be rewritten 
as 


-f I > = < X , ^ I 

+ j s^cz *-2T<v ^ > 


— cO 


(64 j 
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'rt)ls cxpiTCSScs the number o£ particles to be found In as being equal to 

those which would bo there in the obscncc of the potential plus the net num" 
ber (possibly negative) of ortlvols of those streaming in to this interval 
after Jumping In somewhere upstream because of hoving undergone on Inter- 
action. 


In the usual Hilbert space representation of scottorlng, an cquo- 
tlon for the transition matrix may be derived from the equation for the wave 
function by operating on the latter with the potential operator. A similar 
maneuver here yields a similar result; 


((.r) 

In sumtiiary> what we have shown is that the elementary relations of 
scattering theory translate, under the Wigner representation, into elementary 
relations of transport theory. Tlae only thing whch distinguishes quantum 
from classical theory in all of this lies in the details of the Jump pro- 
bability function 


V. CURRENTS, SOURCES AND CROSS-SECTIONS 
FOR MULTIPARTICAL SCATTERING 

Our flrut stop will bo to gcnorall^o the 'density' operator, 
describing o completely locolizcd elementary particle to obtoin a correspond- 
♦ ing operator for bound systems of elementary particles. Consider os a sim- 
plest example a bound state, D , of particles 1 and 2. Lot be the 

Hilbert space ket corresponding to a momentum eigenstate of T) . Then 

where 

42 - C n*. 4?, - ri. 4?^ 1 / c ♦- ) 

Thus is the momentum space wave function describing the Internal 

structure of this system and 

<P I Q> I < 01 D - I (67 ) 

We find it useful to define, as before, 

(ifo) * fdj e (6S’) 

Tliis is now a two body operator but it has much the same properties as before. 

In particular 

^ d p ^ ^ ■= ^ ^ ) 
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1 1 ' ‘ 

I Ro Ko t^’o) I ’ ^ <■‘^ 0 , 'fs>"> , <'■''1 

where fD» -Ro ftttd is the projection operator for \) t Using 

as obvious extension of the expansion theorem written in terms of the element- 
ary particle density operator R'f) which was previously discussed, we have 

Rp i<^o) = J w* '*'<’»• 

• I Ko (ifo) Ri “f.) RiKftJ I Rt-OCi) 

C7Z) 

Carrying through the algebra yields 

where %• * £t > > £.*• the linear combinations of 

appropriate to the center of mass transformation and 

TVp l>, 

is just the Wlgncr representative of the internal structure wave function 
of "0. 

The set of internal phase coordinates of a bound system will be 
given the generic symbol . Thus in the simple case we have been discussing 
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lziry> J2CO .*4>o J ~ 6 A fo 

etc. If the sy8te.> X) were to be InsteacJ | a 3 body system, the set t o 
would denote two internal phase points ond so on. In all cases of 
elementary particles bound into a state we write 

^ p C ^ o) *^pC 

' ' ' ^is ( <^ M ) 

Consider next the integral equation of the stationary density 
describing a beam of particles scattered by a target A in the labor- 
atory frame. Both and 0 may in general be bound states. Assume the 
whole system to be made up of M elementary particles. The incoming 
channel will be denoted by the subscript A 6. The incoming channel 
Hamiltonian is In the absence of any interaction between and (5 

the incoming density operator satisfies the Schroedinger equation 

- -V (. u,.B /’ a* 

or 

Proceeding as before we obtain 

* r r ] (p~t>ao 

( 7 ^) 


26 


for the stationary scattering state. 

To obtain an expression for the cross-section we again construct a 
standard density describing the unperturbed system. 

I { dxo )/■ 

' (7‘l) 

At this point the objection might well be raised that since this operator 
does not really represent a physically realizable quantum mechanical state 
of the system, an error may be made in its use as the unperturbed initial 
state in a quantum scattering problem. A more careful and lengthy deriva- 
tion not reproduced here, shows this not to be the case. It is yet another 
example of a phenomenon noted earlier; the formalism is indifferent to the 
quantum nature of the problem. This density describes a target particle 
localized at the origin with zero momentum while at the same time a steady 
stream of bombarding particles 1^ flows in with momentum and flux 

' 1^6 = 1 1 ^ \ 

Due to the interaction, this distribution is altered to become 

I ft 6 = l ft s > 

.-f f° d* ] CT>-X>ot I (■^oj 

Suppose one were interacted in the inclusive cross-section, A ® ^ ^ ^ 

anything, where C is some, in general, bound state of the system. Asymptoti- 
cally far from the origin, the components of C , if it is a bound state, will 
have negllgable probability of being close together unless in fact they arc 
in their bound state. Thus, at large distances where 

I ^ i Kc > 

is the prob ,lity density for finding C. within d t^c . Note that this inner 
product, ,ce it is in fact the Hilbert space trace of ‘KeCtfe) with the 
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stationacy density descrlbine the collision, includes an integration over 
unobserved final state particles and an integration over ?c.* internal 
coordinates of C . The current of particles C far from the region of 
interaction is 

< '^fcl A (*?2.) 

Proceeding exactly as in potential scattering yields, now for the inclusive 
production rate 

t ^ i Sc-ab ) (^ 3 ^) 

where 

Sc. * Au 

= 'iTc. - Vd x<. < tf(. I A 6 

Note next that, if He is the Hilbert space Hamiltonian acting or the sub- 
space of pratlcles in C. such that 


where is the Hilbert space eigenket of momentum :^c* then it is 

straightforward to prove that 


-i ( He Rc H e") 

= zr<^ ■ 'Kc. iM’. > 

This implies 


Dc ' .IT' - 1 '■<’<•> 





= I R. f 2L- ^ 




(S*7 ) 
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Further note that 

“ o 

where Pc is the generator for all particles not Included In the cluster C . 

This is true because the symbol <^pc\ as used here contains an Implicit 
integration over all the phase coo''dinatcs of tlxcse unobserved particles; 
this fact, in conjunction with equation 12, yields the result. 

Combining equations 87, 88 and 83 yields 

ScjAO -Pc ( fi 

= <:ifa n -t>< -X>r t 

or 

S c. • |!| ft ^ ^ ^ I A 6 ^ (^8 \ 

whBT.’e is the sum of all interactions between members of C and members 

of C. — the set of all unobserved particles in the final state. 

It might be noted that this equation is very similar to its analog 
in the ordinary representation of scattering theory except that it is simpler 
having eliminated all reference to unobserved final state channels by means 
of the closure idea. 

Cross-sections for two or more particles in the final state may be 
derived in a similar manner. For example, consider the cross-section for 
^ ♦-'O anything. The pertinent ti^o particle phase space 

probability density is • Consider surfaces and 

, exterior to and enclosing the region of interaction. The probability 
that a be found in about , and, at the same time, 

that D be found anywhere in the volume V/o surrounded by with momentum 

in ■=*^0 about is 
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Ttie total number of such fairs of particles being produced per unit time is 

oiojAB ‘^•^p 

= ,2!5 

' (air)'-^ <r (()^ ifp I A <''5r' ^ 

Sip^dn 

This can be converted by use of Gauss ' s theorem to 

tTcP-AO ‘ ■*^8 -■ 

=r „£^ jL^ ( <iyj^ taw)'^ 2f<. • ^21*^ 

Vt.— > oo \/o-^ °** 

' 1 dicp Urr)-^ <T if < ‘fe I P 

V/p 

The order in which the limits are to be taken when the integrand is written 
in this form is important. In fact, were the order to be reversed, the 
integral would vanish since for fixed Xp , ^ q>,> i R must 

rapidly vanish as /Vcl —*f oo , and vice versa, due to the necessarily cor- 

.'%r^ 

related nature of the particles’ emission. We can take advantage of this 
fact and add to the above integrand the perfect differential 

TTp . ^/^Xp < L^o\ A 

without changing the value of the integral. This allows us to write 

X « ^ <r^f< 4^0 ^ ^ 
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We have deleted the symbols 


because in this form the order 

( m ' ^ 

= '-(o \ ( 

■< ‘Cc M*.! T>c5 


is clearly Immaterial* Further, 

• ^A2,o ^ < ‘■f' ‘^pB I 


using 


C‘i 4 ) 

or ) 


where C O is the set of particles neither in C nor Q, one obtains, in 
complete analogy to the one particle inclusive case, 

^P;/IO ^ J ) *^0 ^ 

- ^ Sco-A^o ^ J 

where 

Sco;(V0 ‘fo> = < <^C M’dI T/o \ A B“^*> a?3 

X/n is the sum of all interaction between members of the set CO and the 
unobserved set CV as well as interactions between members of C. and members 

of P . 

Thus we have obtained the general Integral equation for 
and shown how to x^rlte down the expression for any desired cross-section. 

The formulae are simple generalizations of those obtained in potential scat- 
tering and analogous to those of ordinary scattering theory. 


J 



VI. MULTIPLE SCATTERING EQUATIONS 


The analog of the Lippraon Scliwinger equation for bhe multiparticle 


problem, given in equation 80, may be rewritten as 

“ |fVB> CgAt^^TAD \ ^ 

where ^ 

^ di:e“ ^ 

« — 

= - C Dab -7 ) ■' ■ -I ^ *0 


as ) 




- -P - V^B ooo) 

are the Green's function and interaction operators respectively. As in the 
case of the Hilbert space representation of the problem, there is an alter-* 


native form for this equation: 

1 A 6*^^ > - 1 A B > -ft- C) Xa3 \ A (#o ' ) 

where 

C 7 - - Cx>- 7 )"' 


The full Green's function Of may be expanded in the usual manner. We write 
the interaction as a sum over pairwise parts 


x> - 


(lot.) 

X. " 

£T 

('/OJ) 


M 
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where << runs over all p.nlrs of particles in the system Now define 

a 'jump' operator by 


where 


3*e< ~ O «5* 


C) o - - C "Oo y ) 


0 

( 1 o jr \ 


Tlius the two body jump operators take the place of the two body transition 
operators in the usual development. There is one interesting difference 
between the two representations however; neither the Green's function or 
consequently the jump operators used here have the energy of the ^ B 
system as a parameter as is the case with the transition matrices in the 
Watson expansion. Tlius there is no possibility of an 'off-shell' or 'half 
off-shell' type jump operator and, in fact, the same jump operator is valid 
for any energy and any system in which the corresponding two particles appear. 
This fact does not effect the algebraic structure of the equations which Is 
the same here as in the usual analysis so that we get finally 


Q - ^ C-Jp X Ca 

T - X t: 


6,0 X , 


cr ^ - X. 


2. 

J _ 


T G,c X. 


Xu + 


S- X Cj. T. 


(. f O t ) 
do-,) 
( loe) 

CIO*? ) 


Q - (sjo Xei +* 


±pc *■ 


6^0 6,0 

Cio L T'*'- x^Ti 


T‘-> 


(HI ) 
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Tlic last result Inserted Into equation (101) yields 

( A 1 A »5 > ^ ^ G)o I A 

where the prime on the summation Indicotcs that o< la restricted to all pairs 
not InteractinB In the AiS channel. The perturbation expansion o$ equation 
(110) inserted into .equation (112) is the multiple scattering equation? 

\ P ’ I R o •> 

•- 2.' c-to I 

<4 

^ ^ Qo T« Q„T.< \RB> 


Cii^) 


( »o 3 


VII. PROPERTIES OP THE JUMP OPERATOR 

Before extending further the multiple scattering theory begun in the 
previous section, it will be useful to examine the properties of the elementary 
operators appearing in the expansion given in equation 113. 

Consider first po tential scattering and the matrix element 

< If * I i (f > “ < I f ^ y 

r •=** 

- J ^ f c M4-) 

o 

According to the probability interpretation of tlie symbols developed in tltis 
paper, C I is the probability density of a particle to be at 

If ^ at a. time 't, after it was localized at ff » Tlie integrand in 114 is thus 
the divergence of the current set up at h at Applying the equation of 

continuity we get 

ftr ! Cfy *: L ^ } 

- J ir SUf^^ir) ■*- ^ (p > 

where S^v>) describes the sources due to the interaction and 
the particle probability density at any time <5r. Since is a time 

rate of production, its time integral equals the net production. We can now 
describe the situation as follows: An external source pumps a net of 1 

particle into the system at ^ . The system over a period of time reacts and 
redistributes this source distribution from ^ V to • 


34 


35 


Since CJ“ is in the nature of a probability one must hove 

This is in fact satisfied by the theory already since , using equation 106, 


<T ^ \ V T 6) 


and 


(I / t ) 
O / V ) 


I d tf ' < I J « o 

Identically. 

The relation between the operator C” and the function already 

Introduced in the discussion of time independent potential scattering (equation 
60) is, 

Where ^ L b^- ^ . 

b|? is the impact parameter two-sector perpendicular to the entering momentum 
, and 2:^“**^ “ along the beam line, where the accelerator is situated. 

The equation says that the accelerator is a source distributed randomly over 
the plane •=? — oo with an average density of one per unit area. Tlic sys- 
tem reacts to this 'external probe, by redistributing this density. Thus a 
particle introduced at f<r will reappear at with probability density 
I tr f • Tlie cross-section is 

0-(^, fV') = a ‘fZ 

In principle o“ could be computed using 
< <f’| CT| 

_ ; H -t 12 . . , V ^ I 


(11? ) 


I e 






if the spectral decomposition of H is Icnown, This seems to be a difficult 
task. It will be seen that one can often avoid such an explicit computation. 
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Tlie classical limit Is of practical Interest and can bo discussed 
In some generality. We have, as discussed In. Section III, 

X) ^ 

< £('' 1 XTc \ Cp =• Jtr ‘f 7 e I 

= I if 7 


To see what this says lmagii\o a numerical, integration of the classical equa- 
tlo*.! of motion. Time is segmented; the first Interval is O 't' < -i*/ ; 

the next is "tt ^ t: < ^ etc. During the "i ^ interval the trajectory 

is approximated by a straight line .* 

%{ iti-t = 2.-.. '"I 


while the momentum is constant. The momentum changes d is continuously, 

as if by an Impulsive force, at the Instants * • H^us 

de £/> 7 e /'/’,• > 

*1-1 

- i/y 

t C'- <v7?/«, #’.• > •- ^V7^,v,, :f,v, V] 

Ctzj ) 


having used, 


T>o ir 




C, ' T ^ i ! 


‘Zf ,-Pf > 7 * :t y T-. t’i 
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Thus the reaction of the system in the classical limit to on imposed 
unit source at is a line distribution^ along the classical trajecto>ryi of 
sources and sinks. At each spatial point subsequent to X there is a unit 
sink for the old momentum going into the point and a unit source for the new 
momentum leaving the point. 

Another limit of considerable practical and conceptual Interest which 
can be discussed in some generality may be termed 'dilute'. This will be 
applied to a multiple scattering problem when the distance between successive 
collisions is large compared to other oiaracteristic lengths of the problem — 
wavelengths and ranges of potentials of particles. In particular this limit 
Implies no necessary ordering beb'^een wavelengths and ranges whereas the 
clasaleal limit requires wavelengths to be smaller than all other characteristic 
lengths . 

Consider, for simplicity, scattering by a system of fixed scatfcerers 
centered at points The multiple scattering series is still given by 

equation 113; the subscripts now range over fixed scatterers* indices. In 
any one of the terms of the series in which the <=( ^ jump operator 
appears it operates on some particle density function; call it i/>> * Xn 
the dilute limit, Is set up by sources centered far away from . 

This leads to the conjecture that any such density must be slowly varying over 
the scale of distance set by the size of the ~ scotterer, the range of 
its potential. We could then, expand about ^ h and, in the most favor- 

able case, keep only the lowest order term; thus 

J <4 tf T«< I 4 ^ ^ I ^ Cl Z-4 J 
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Howovor> Quantum Machanlco Implies the existence, In principle, of 
substontlal oscillations In associated with interference* These would 

vary un the ncale of a imvelen^Ui and contradict the ccm.loot are. Thin ipu-ition 
in rC'iolved in the furtlier quantum miT.luinij'aj ri qu'n ;”cn\t tliaf. nouret-n !-e iiln" 
tributc<l ovi-T a certain min imum i’c 'ion in phanc: ii,iaco. Oar i1 Jatlonn <hin to 
-'Untrlhuted ;;ourccn tmnh out Increanlnf'ly v;ll.lj diJit.nue.-in the dilute Kt; i t 
they dloappenr. 

Wlu'iv tlie condition of d’lutfneuo !■; rtnt ntronn ciViiipji to all.fy all 
worriun annoclated with the approximat ..on «f i-quation 124, other factorn may 
nnrjt'it. RandontnenG In often Invuhetl to aecoiupllrdi phase averaninn. A rela- 
tively fit-uKiity packed medium nuch rr. a liquid or a nuelt-us will ponnetui only 
chort r<anfic order. The approximation In equaidon 1.24 I'r'rht he >yjod provided 
only that contaiiu; no contribution in vdtlch a nt .ir te.’fphbor of wau 

th‘‘ lar<t ncattcrer. Temfi in the iiuiltiple rcattf>rlnu serier. whic!. doncHbet 
fatacesfilva acattcrlnRa by nrar nelj’hbcjr.o mip.ht be treaii'd si-parately, leadini’, 
to an expansion in powers of correlation fimcfcl>>juj of the i t lUuw. 

Another escape hatch opens If the uav vleni'th *'i (be jinrtJ<’l('-5(!:) 
belnp ficattc'red Is short compar(.‘d to a lon^ ranj;o component- of tfie potential:*, 
doi.nn the scattorinj;. It Is stralpht fon;ard to show that the multiple saMt- 
terJnH series can be rc~wrltten so that the elas.slcal pnipap.ator C^c replaces 
^4o « Concurrently the Jump operator Is replaced by one wiih h represent:. Ho 
stochastic quantum jump correction to tin classical motion .as discussed in 
Section III. This correction is due mainly to the shortraiiRC component of the 
potential. .Such cores ray be spaced v/ldely enough apart for the dllutencss 
condition to apply. 




j 
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terer 


is the source distribution created by the scat- 

when >r impinges. This must always be integrated over another 
density, say calculate an obiservable probability, 

But the arguments applied to , apply to , using time 

reversal invariance. Thus one can write 

5 dtf 2i'< ^ oar \ 

Combining equations 124 and ‘’.25, 

I if> = SC2i'-Z-) S’tz-z,-' 

• i < ‘f’l O’-*-) 


In contrast with the general expression for 0” given in equation 
120, this approximation for the closely related jump operator is easy to 
evaluate and interpert. Note that 

is the density for a unit concerning plane wave state. It can therefore be 
written as 

C 2 IT I ( \ ^ ) 

where ip> is a Hilbert space ket. Next note that 

Tm = c? X- ’ 

which, in combination with our integral equation for scattering, equation 101 
yields 

= (IIT)* t),'' Q I 
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Hence, 

= fciir)^ j <• T< ' <2;:g'l G»o"'[,l ( 

- Ui f ><f 

= (^ir>^ j dTj' y. a/o2s' [;<■ ^‘*’l Rtz'^') I 

- <-f 1 Rti', ^'>\ 

using equations 5, 7» 18, and 19. Inserting an expansion for 
such as given in equation 9 and writing 

< I 

+ ( fe<.^)+.o X 

< ip'*"' 

“ c -^L I dP - i ? 

^- Cfe(^)-iO - fe( ^ 


Clio) 


Ci3i ) 


(isz) 
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for the usual expression relating outgoing wave solutions to Green's function 
and transition matrices, one gets, 

i dl' dr Tv I T, 

^ errr)^ j iiz' e'^'^ 

' / )“’X. 

+ ( 1 ) 4-10 - fc <^' 4 -i^'U"‘(fc -7 o - fcf ^ '-iq.') )'* 


The terms linear In / are the complex conjugates of each other. The first 
of these Integrates over to become 

i%TTy^ jdr' e ' ipt z’- s') C(:tsy*'‘^ - 

C/3+J 

with — \ ' The integral over c(r' will give a J~-function 

distribution in so that we can write inside the integral 

- (rijp *-q') s. -q'. ^ - s' T ■ 

Integration over /J'jc* gives 
The sura of ^4ae linear terms in T is 

^ U <■ S ; ^ ) . 


-X l3.n) ^ Sc^'-s ) 


0J7) 
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The quadratic term may be easily evaluated by noting again that the ^ ^ 
integrotion restricts the Integrand to the neighborhood « a , Thus we get 

It, C:^;£fr 

' - e ) ~* C -»'0 - ^ )"' 

^ t-Z.Tr\^ fciTT) \ ^ ( 1 3 >^ ) 

Hence we obtain in the dilute limit, 

< cf W / v> > 

( 13*7 ) 

Using the optical theorem and the relation between the transition matrix and 
the differential cross~section 

V" t-p) ^ --Z. tiTTl"’ T(^/ 

'tr \T(:^'i^)r cm) 

where related to the differential cross-section by equation 

61 , this result can be rewritten as 

< ii>'i X< / 

= S(TC''--Xrf ) S’fTC -2-*) 

^- -iT CTL (:p'> j 


(I 4 Z.) 
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This equation says something which could have been written down 
right away using elementary considerations. In particular, a density 
Impinging upon the fixed scatterer gives rise to a source distribu- 

tion function 

(T «r)* ) 

043 ^ 

The first bracket is the incoming flux per unit area of particles between 
^ and at . This is multiplied by two factors in the second 

bracket; the first gives the number produced per unit ^at posi- 

tive source term; the second gives the number lost to the incoming beam— a 
negative sink term. 

This result can be Immediately generalized to Include non-fixed 
scatterers. Let d now refer to a pair of particles, e.g. particles 1 and 2. 
Then we have 

ifi' \ l 

-S sc X.' -X-') "SC P..'-P.x) Sc^.ly yc?..) 

C7.- OH) 
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where Xtt.> refer to phase coordinates of the center of mass while 

apply in the center of mass. This result may be obtained from the for- 
malism but can be written down directly on physical grounds. 

The advantage gained in Iteving derived this approximation via this 
formalism js th«t it is now part o^ quantum mechaiaical scattering theory. The 
possibility of systematic refinements to any approximation is tliereby created. 
We have already discussed quantum corrections to the classical limit in Sec- 
tion III. The possibility of combining the classical approximation (for mo- 
tion through the long range part of a potential) with the dilute approximation 
(for transitions induced by short range potential cores) has been mentioned. 

In addition, finite size corrections to the dilute approximation, arising 

when the range of the potential is not negligably small compared to the dis- 

[81 

tance between scattcrers, have been calculated.^ ^ The relation between the 

dilute approximation and Monte Carlo calculations has been discussed else- 

[ 5 ] 


where. 


VIII. MULTIPLE SCATTERING EXPANSIONS OF 


INCLUSIVE CROSS-SECTIONS 

The previous section developed certain limiting cases which are 
useful for future applications and serve to Illustrate and clarify the mean- 
ing of the basic Jump operators that appear in the multiple scattering expan- 
sion. We now return to more general developments of the theory. Formulae 
for the simplest inclusive production processes are rev;ritten in a multiple 
scattering form which facilitates their interpertations and computation in 
many cases of interest. 

It is useful to define a symbol, ' J , for the sum of all 
scatterings beginning with the o< and ending with the p pair. Simple 
expansion in powers of the Jump operators will verify that it satisfies, 

X- ^ ^ O - Cjo 

Furthermore, using equation 106, 107 and some strictly algebraic manipula- 

« 

tions, one can show that 

T = Tu * Tfi G, Xw . 

Particle ^ impinges on a composite target A. Consider first the 
Inclusive source distribution function for B observed in the final state. 
Application of equation 89 shows. 
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^ I -^B ^ AG*'**'> 

= JXo c I + . (HtJ 

Since I© - 3!^ 5 , the sum of alt interactions between and constituents 
of A » ons Iminediately obtains 

So.,m < '^o' I *■ ®“-'’ T I A B > ('Hr) 

where a and Ca.' run over all constituents of If, for example, A were a 

two body bound state composed of particles labelled 1 and 2, the multiple 
scattering expansion of equation 148 yields 

+ '3 01 Crjo 0x 

These expressions will be evaluated in the dilute limit in the next section. 

Expressions for the inclusive production of a single particle other 
than the incident one (knockout) , follow in a similar manner. These source 
functions are, of course, to be Integrated over space to give the total pro- 
duction rate or, equivalently, the cross-section. 

Several examples of practical Interest Involve two or more particles 
observed in the final state. They may be in bound or unbound condition. Our 
intuition would like to relate, for example, a two particle Inclusive produc- 
tion rate to that of a bound state of the same two particles. The present 
formalism brings this out in an especially graphic manner. 
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The Inclusive source function for the bound state production of 
made up of constituents 1 and 2 originally In A is given by equation 89: 

^ (iro; 

that for 1 and 2 unbound Is given by equation 97: 

St tf an W/ ) ■= < ‘V/ tf/i :r/r. I PS Ci ^ ' J 

Mov 

x» = 1,0 » ^ x,«,. 

a^t,iL 

+- Tx G ^ Xxc c»^2 ) 

while 

x.l. = To + X,,. C/j-3) 

Thus we see that the basic source functions required to calculate bound and 
unbound production differ only in the term 

^ I X IT, \ AC ^ 

This terra represents final state interaction between the observed constituents. 

It is the sura of all contributions to the multiple scattering series in which 
the very last interaction is betoeen the constituents of the observed final 
state. This final state interaction must be deleted in the computation of 
bound state production — incJ.uded in the unbound. 

Except for kinematics which especially emphasize the role of the 
final state interaction (e.g. at low energies and low relative momenta 

- •i^t. of the observed pair), this term is unimportant. It is useful to 
define a source function which neglects this term] 


S,i-Ao X„' I 
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Then _ 

- topt-jc', ip'; §■ .-..j ^ „ ( cf,' iK' ) 

S IJL M O ^ ) 

where ^x^A*5 cross-section less final state interaction contri- 

bution. In the above formulae we have used: 

A ”X/ Jtr = dK ‘ Jyf,' 

^o' - -P.' * , 

( Mr t(j/ - I'll Jpr') /(M. VMr) . ^ 

These equations say that the source function S*, giving the spat- 
ial probability distribution of the origin of the particles produced in the 
interaction already contains all the Information necessary to compute the 
bound state production rate. One just convolutes S with the Wigner density 
describing the bound state. Thus, for example, a Monte-Carlo calculation 
approximating the multiple scattering series and hence providing an estimate 
of S , can, without further effort, provide an approximation to the bound 
state production rate. 

If one further approximation, often valid in practice, is made, 
these relations become even more striking. The dependence o£ S on is 
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assumed to bo much less rapid than that oC 7i/|j which is peaked ot . 

In this case one has 

t‘4?o) lie « ^ V | M'd 

TU I ^ "Tl ' 

)'^0 * J (art)'’’ <) Tc/ i3ir)'^ 

^ 11.^ Al o / -t: ; 2t‘ / -TT > 

The Integrands differ only in the substitution of the spatial probability 

l'M'o«i'>r= wof^'^g') 

in one for a factor of other. Thus the ratio 

Op|A& CJ^O ) / CT,2./».a ( % y C*" ) 

is related to the volume from which the particles 1 and 2 emanate and hence 
the mechanism of their production. 
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IX. MULTIPLE SCATTERING BY A TWO BODY 


BOUND STATE tn THE DILUTE LIMIT 

The multiple scattering expression for inclusive scattering of B 
by a two body bound state A is now evaluated in the dilute limit. The 
methods and results of this formalism mny be connected with others in the 
context of this well known example. 

The contribution first order in ^o»t to this process is given 
by (see equaiton l49), 

^ I Ttj, ( A o > 

' < <fa I Xei \ > . c 

Since we are in the laboratory system •= O . As before 

are the relative coordinates and momenta of 1 and 2. 

The dilute expression for can be written in the present 

notation, using equation 144, as 

<: i^e,' { Se, I 4>, cf^y 

- c fei I ioi'y Scrl> 

[ •- -fe.) 

* ■e.', 


( ) 
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5 61 " ^ 

arc the phase coordinates of the 6^ center of mass; 

<^f,\ - o • 

refer to relative coordinates » specifically 

%Qt ~ "K 6 - y I 
and 

<^7, c^op, which really only depends on -jPg , , Is the total cross 

section for ^ on 1; C^, . which really on depends on 

/ A ^ 

anti > Is the differential cross-section for scattering into 0^e \ • 

Its relation to the usual center of mass differential cross-section is given 

by equation 61. Combining these equations and integrating over spatial 

functions yields 

J d dif Last) < 3S, ^ ^ 


52 


where now 

^ 

rhe first term gives the contribution to the total croos-scction due 
to collisions of 6 with 1 only. If we add to this the corresponding contri- 
bution from 2, the firut order estimate of the total cross-section is t 

0 

Tlie second term in equation 163 is the well Icno^vn single scattering term 
originally obtained by Goldbergcr, 

to which one should add the corresponding term for scattering of on 2, 

The contributions second order in the jump operators provide shadow 
corrections to the first order results as well as double scattering. The 
contribution from the ‘1 0 « term in equation 149 is 

<r I Oftx t?o Toi 1 4*. ifi - (ii-^ ) 
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Expanding over intennedlato states and using 

where 0‘-7 ) 

Ifo (<f, <r* = 2 *ir* , ^ 

we get I 

< I 

■ < ^," >-e." tfo I To./cf. If. > 

- M”* 1 Tb, I tf, CfB> 

(ity) 

Inserting this result into equation 166 and integrating over the spatial 
functions, we get 

C dt f o (i!n'3-’Wrtfz,.e> 

• ‘Z ' r*- * + iTo" ’ f ^ 

Sc^e'.: 

■ I ;% -5, 1 S( ^ffV - iPe -^ ) 

■ [- i>) 5 -:po.3 ^-<3751 ^/> 1 

1 

(I ‘O 3 

i 

* ‘1 

i 

■\ i 
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The spatial S function requires that at some time t after collision of i3 
with 1, the position of B , "Xj , coincide with the position of 2, 

'JC.x lT . The only new notation used In equation 169 Is for 

the relative momentum of 1^0 and f, , -f «x for that of and , 

for that of and and, 4?o i for that of fe and . 

The term quadratic In CT represents a shadow correction to the 
total cross-section; the terms linear in Cl" represent shadow corrections 
to single scattering; the term quadratic In the differential cross-sections 
represent double scattering. We consider each In turn. 

The shadow correction to the total cross-section given by this 
expression Is 

- ) cxwr’ £ dt i^h *2TI L-lVfi ) 
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which must bo added to a similar term coming from collisions in the reverse 
order to get the total shadow correction to second order. This expression 
simplifies somewhat in the high energy limit where it can be compared to the 
corresponding correction due to Glauber. Using 

\ TTg ±iy { ~ -iTe 
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for the bound state spatial density, we get, 

Taking into account that the contribution to the shadow correction 

coming from the reverse order of collision, one finds that the total shadow 

171 

correction in the dilute limit is tv ice that obtained by Glauber in his 
simple approximation to black sphere scattering. Tlie reason for the differ- 
ence in these results may be traced to the fact that they each hold true in 
different, non-overlapping, regimes. Our result is accurate in the dilute 
limit in which the Glauber approximation is invalid. It should be re- 
emphasized at this point that the basic fonnalism — the Wigner representation 
of scattering — is not limited to the dilute limit but is susceptible to 
other approximations which may lead to simple formulas valid in other 

regimes . 

The shadov; correction to single scattering is composed of two parts. 
The particle casting the shadow may li» either between the accelerator or 
between the detector and the single scatterer. In either case one gets a 
reduction in the single scattering — a ’unitarity correction' since it is 
merely a, manifestation of probability conservation. The shadow correction 
to single scattering by 1, to be added to the uncorrected term given in 
equation 164 is obtained by straightforward application of the equations 
just discussed yielding again in the dilute limit, 

- ^ dt: U'n')"'* 

I 'T/i - -IT I / , 

. / / (Te, 

+ I T»'-£l Ho 

-^3 


{ /7T- ) 
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where is again given as in equation 164. A similar term reduces single 

scattering on the other particle, The physical origin of these terms is 
transparent. 

We might note that the effect of sucli corrections is quite striking 
in reducing the backward peak observed in proton-deutron elastic scattering 
( 50% correction) . It may be that a simple formula such as derived here 
(perhaps somewhat refined to take into account the not completely dilute na- 
ture of nuclear matter) will be quite accurate in calculating such shadow 
corrections for many nuclear targets, including the deutoron. This would be 
especially helpful at those intermediate energies at which the Glauber approx- 
imation begins to break down. We also note that, as before, this formula 
simplifies considerable when Vjj is much greater than 'V ~ , the Fermi velocity. 

Finally, the double scattering term obtained from equation 169 is 

\ (aiJ")-’ dsi' £ Je -vi" 

. I fa -f/ (•fe'tF.'i C/7JJ 


where now momentum conservation gives 

- -fa -Jp - 

To this must be added the corresponding term from the reversed order of col- 
lision to get the total double scattering. 
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We have here a 5-dlmenaional Integral since there are two energy 
conserving 5^- functions implied in the clastic cross-sections. Wlien 'iFg 
is large compared to t/" , we can again effect a considerable simplification 
by integrating analytically over One is then left with a 2-dimensional 

integral. In cither case this integral is especially suited to Moute-Carlo 
integration techniques. 
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FIGURE CAPTIONS 


Figure 1: This illustrates the first order correction to the classical pro- » 

pogator given in Equation 41. A particle is known to be at ^ 
initially and one asks for the probability of it being at if ^at 
a time later. In addition to the classical value of this pro- 
bability there is a contribution corresponding to the particle 
traveling a time - fc ^ along its initial classical trajectory 
and then performing a quantum jump to a new classical trajectory 
which brings it to ^ in the required time . The trajectories 
are pictured in position space emphasizing the fact that, for a 
local potential, there is a jump in momentum only as with a classi- 
cal stochastic impulse. The probability that the particle makes 
the jump is given by Equation 41. 

Figure 2; This illustrates the classical limit of scattering in terms of 

source functions as given in Equation 123 and discussion subsequ- 
ent to it. A circled portion of the spatial trajectory is 

enlarged. The continuous force is approximated by a sequence of 
inpulses. A particle starts at with momentum . Thus at 
7^1 thexa is a unit source of such particles. At 7^ receives 
an impulse, changing its momentum from to Thus at Xu 

there is a unit sink for and source for , and so on. In 
this picture, classical motion is a consequence of the potential's 
creation of a line of sources and sinks in response to the extern- 
ally imposed source at ^ J^, . 
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Figure 3: This Illustrates two phenomena generated by the second order term 

o£ the multiple scattering expansion for scattering of projectile 
by a two particle bound state A (See Equation 149). When 
hits constituent ^ it happens, in this instance, to have momentum 
and be displaced from constituent 1. by_2^ . The solid line shows 
the trajectory of B . Between collisions B , having gotten inter-* 
mediate momentum , dldplacos ItselC along while 

constituent 2. , since A was inltally at rest in the laboratory 
frame, has momentum and has moved -jy”"*. The collision be- 
tween B and i at a time "t after the first collision has two 
consequences; First it reduces the number of (single scattering) 
events into Scf » i.c* particle 2. shadow events shining from 
particle * . This phenomenon is described by the first bracltetcd 
term in Equation 172. Second, it increases the number of (double 
scattering) events into . Tills is described by Equation 173. 





